Power-law tailed spectra from equilibrium 
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We propose that power-law tailed hadron spectra may partially stem from a matter 
in an unconventional equilibrium state typical for non-extensive thermodynamics. A 
non-extensive Boltzmann equation, which is able to form such spectra as a stationary 
solution, is utilized as a rough model of quark matter hadronization. Basic ideas about a 
non-extensive simulation of the QCD equation of state on the lattice are presented. 

1. Particle spectra 

Statistical models have been often applied to hadron physics. Starting with Rolf Hage- 
dorn's statistical model of meson resonances [ Q], several attempts occurred to describe 
hadron multiplicities in elementary collisions by means of statistical distributions. The 
very idea of a phase transition between confined and deconfined quark matter relies on 
traditional equilibrium thermodynamics. The search for quark matter began assuming a 
local thermal equilibrium in an otherwise exploding fireball [13 El- Experimental particle 
spectra are, however, not purely exponential: both exponential and power-law regions 
have been observed in pion, kaon and antiproton spectra. Although the traditional ap- 
proach explains the power-law tail at very high px values by pQCD calculationsf 4J, the 
non-extensive statistics provides a unified view for the whole spectrum. Pion spectra from 
heavy ion collisions at RHIC seem to contain a power- law part exceeding the scaled pQCD 
yield Another point can be made by inspecting the minimum bias pion p^-spectrum 
form RHIC AuAu collisions at 200 GeV (Fig.l in Ref. A non-exponential fit can 

already be made at the pr-region between 1 and 4 GeV. The extrapolation of this fit 
almost coincides with the fit to the whole observed range between 1 and 12 GeV. So one 
concludes that the power-law behavior is not only a very hard scale physics. Furthermore 
in the non-extensive statistical approach there is a connection between the soft properties 
(temperature T) and the hard ones. 

The experimentally measured specific hadron spectra may reflect statistical properties 
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of the precursor matter [Uj. Fortunately transverse momentum spectra are influenced only 
partially, at their low end, by final state interactions and late resonance decay[|H]. Since 
the relativistic energy is given by E = mxcoshy with transverse mass mi = \jvr + 777,2 
and rapidity y for a particle with mass m, the best way to study statistical equilibrium 
distribution of hadrons is the comparison of m^-spectra at rapidity y = for different 
particles. A universal behavior [Hj indicates that the one-particle distributions depend on 
the energy only and not on all momentum components: a basic feature of generalized and 
conventional thermal distributions. 

2. Non- Extensive Boltzmann Equation 

Non-conventional distributions can be based on a non-conventional entropy formula, 
which replaces the Boltzmann entropy. Such a formula is the Tsallis entropy, discussed 
vividly in recent years. This non-extensive thermodynamics is intended to be an ef- 
fective theory for non-equilibrium and long-range order phenomena [ E3- Its canonical 
distribution is a power law, which occurs in particle and heavy-ion physics experiments. 
As nonlinear models, two generalizations of the Boltzmann equation have been investi- 
gated: The generalization of the product rule for probabilities (dropping statistical inde- 
pendency) leads to a non-linear Boltzmann equation [ llj, while considering two-particle 
energies composed by an extended addition rule mounds in the non-extensive Boltzmann 
equation [ 12. 

The general structure of the Boltzmann equation describes the evolution of the prob- 
ability fi — /(pi) of a one-particle state by considering possible transitions to and from 
other states: fi = J Wi 234 (^34,12 — /i2,34)- Here the dot denotes a total time deriva- 
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tive (Vlasov operator) comprising the essential evolution of the one-particle phase space 
density, f\. The indices 1234 refer to two particles before and after a microcollision. The 
transition probability, W1234 contains conditions on conserving momentum and energy: 

^1234 = M% 2U 5({pi + p 2 ) - (p 3 + P4)) S(E 12 - £34), (1) 

with E12 total two-particle energy before and £34 after the collision. The particle density 
factors, /i2,34 and ^34,12 weight the transition yields for a 3 + 4 — > 1 + 2 and for a 1 + 2 — > 
3 + 4 process, respectively. In our approach[ 12] we keep the statistical independency, 
/i2,34 = /1/2; but generalize the energy addition formula to a nontrivial composition 
rule E12 = h{Ei 1 E 2 ). Rules not being a simple sum, h(x,y) 7^ x + y, present a non- 
extensive energy composition. Associative rules can be mapped to the simple addition [ 
H3j: X{h) = X(x) +X(y), unique up to a constant factor. In each microcollision X{E\) + 
X(E 2 ) = X(E 3 ) + X(_E 4 ) holds, and the stationary solution is hence given by f(p) = 

exp(— X(E)/T). A statistical dispersion relation is obtained due to the mapping of the 
general composition: the total sum, X(E tot ) = J2iX(Ei) is conserved. 

In order to utilize such a non-extensive Boltzmann equation for the stationary state 
Tsallis-distributed parton matter, we consider an evolution from initially boosted Fermi 
spheres. During the simulation we take out pairs of partons with a color singlet against 
octet probability 1/9. This rate turns out to be low enough not to disturb the Tsallis 
distribution remarkably. Fig^ shows the resulting energy (tut) distribution of such pairs. 
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Final hadron distribution (N col | / N = 4.5, m = 0.14 GeV, p = 1 / 9) 
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Figure 1. Pion spectrum emerging by sampling with a probability 1/9 from the NEBE 
parton cascade with massless partons. 



It is interesting to note that an exponentially growing mass spectrum, originally pro- 
posed by Hagedorn and recently checked again latest experimental data in Ref. [E], 
with its famous consequence of having a limiting (or Hagedorn-) temperature for such 
a system, can be reconstructed on the basis of Tsallis distributed quark constituents. 
This approach [ITo] assumes that the Tsallis distribution of the quarks and antiquarks is 
folded into mesonic and baryonic distributions of the conserved total energy satisfying 
X(E) = J2iX(Ei). The Hagedorn temperature is given by Th = E c /d. 

3. Lattice strategy 

The implementation of the Tsallis distribution in lattice field theory can be based on 
the superstatistics approach [ [TH]. The Tsallis expectation value of an observable A[U] 
over lattice field configurations U is of interest. It may include the timelike link length, 
say on the power v: A — 8 V A. The Tsallis expectation value then is an average over all 
possible a t link lengths according to a Gamma distribution of a t /a s . We obtain: 

(A) T8 = T^T^T ldee c - l e- c6 I VUA[U}9 v e- s ^ (2) 
^ts 1 yc) J J 

with Z TS obtained by requiring (1)ts — 1- The 9 dependence of the lattice gauge ac- 
tion is known long. Due to the time derivatives the electric(" kinetic") part scales like 
a t a? s / '(a%a%) = a s /a t , and the magnetic ("potential") part like a t a z s j \a 2 s a 2 s ) = a t /a s . This 
leads to the following expression for the general lattice action: S [9, U] = a 9 + b/9, where 
a = S SS [U) sums space-space, and b = S ts [U] time-space oriented plaquettes. In the c — > oo 
limit the scaled Gamma distribution approximates 5(9 — 1), (its width narrows extremely, 
while its integral is normalized to one), and one gets back the traditional lattice action 
S = a + b, and the traditional averages. For finite c, one can exchange the 9 integration 
and the configuration sum (path integral) and obtains exactly the power-law-weighted 
expression: (A) TS = J VU W v>c [U]A [U] / JVUW , C [U], with the Gamma fluctuating 
time-link averaged general weight factor, 

Wv c = [ dee *+ c - l e- c e e- s ^ . (3) 

Tic) J 
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The 9 integration can be carried out analytically using the replacement 9 = e* yb/(a + c). 
The result contains the K Bessel function: 

The K-Bessel function has an exponentially decreasing asymptotics, so we are in principle 
able to utilize known Monte Carlo techniques in order to calculate Tsallis expectation 
values. On the other hand we cannot simply use old data, produced according to the 
weight e~( a+b \ because the argument of the K-Bessel function is not a + b. This makes it 
necessary to redo lattice calculations - but only with a slightly increased effort. 
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